Abelian anomaly and neutral pion production 
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We show that in fuUy-self-consistent treatments of the pion; namely, its static properties and 
elastic and transition form factors, the asymptotic limit of the product Q'^G^*^^o{Q'^), determined 
a priori by the interaction employed, is not exceeded at any finite value of spacelike momentum 
transfer. Furthermore, in such a treatment of a vector-vector contact-interaction one obtains a 
7*7 tt" transition form factor that disagrees markedly with all available data. We explain that 
the contact interaction produces a pion distribution amplitude which is flat and nonvanishing at the 
endpoints. This amplitude characterises a pointlike pion bound-state. Such a state has the hardest 
possible form factors; i.e., form factors which become constant at large momentum transfers and 
hence are in striking disagreement with completed experiments. On the other hand, interactions 
with QCD-like behaviour produce soft pions, a valence-quark distribution amplitude that vanishes 
as ~ (l — x)^ for a; ~ 1, and results that agree with the bulk of existing data. Our analysis supports 
a view that the large-Q^ data obtained by the BaBar Collaboration is not an accurate measure of 
the 7*7 — >■ Tv'^ form factor. 
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I. INTRODUCTION 

The process 7*7 — > tt*^ is fascinating because in order 
to explain the associated transition form factor within the 
Standard Model on the full domain of momentum trans- 
fer, one must combine, using a single internally-consistent 
framework, an explanation of the essentially nonpertur- 
bative Abelian anomaly with the features of perturba- 
tive QCD. The case for attempting this has received a 
significant boost with the publication of data from the 
BaBar Collaboration [1] because, while they agree with 
earlier experiments on their common domain of squared- 
momentum-transfer [2, 3], the BaBar data are unexpect- 
edly far above the prediction of perturbative QCD at 
larger values of Q^- 

Herein we contribute toward understanding the dis- 
crepancy by analysing this process using the Dyson- 
Schwinger equations (DSEs) [4-9], which are known to 
have the capacity to connect nonperturbative and pertur- 
bative phenomena in QCD. In particular, the connection 
between dynamical chiral symmetry breaking (DCSB) 
and the Abelian [10-16] and no n- Abelian [17] anoma- 
lies is understood, as is the manner through which the 
perturbative QCD results for the large-Q^ behaviour of 
transition form factor can be obtained [18, 19]. 

As part of this analysis, we will elucidate the sensitiv- 
ity of the 7*7 7r° transition form factor, G^»-y^((5^), 
to the pointwise behaviour of the interaction between 
quarks. We will use existing DSE calculations [20] of 
this and the kindred 7*7* — 7r° form factor to char- 
acterise the Q^- dependence of Gj*^Tr{Q'^) which is pro- 
duced by a quark-quark interaction that is mediated by 
massless vector-bosons. For comparison, we will compute 



the behaviour obtained if quarks interact instead through 
a contact interaction. Such comparisons are important 
to achieving a goal of charting the long-range behaviour 
of the strong- interaction in the Standard Model [21]. 

In Sec. II we describe a symmetry-preserving regular- 
isation and DSE-formulation of the contact interaction, 
and explain how a dressed-quark comes simultaneously 
to have a nonzero charge radius and a hard form factor. 
In Sec. Ill we discuss the 7*7 — )■ tt" transition form fac- 
tor in detail, describing: its connection with the Abelian 
anomaly; its asymptotic behaviour in QCD cf. that pro- 
duced by a contact interaction; and how the nature of the 
interaction determines the pion's distribution amplitude. 
Section IV places our results in context with extant data; 
and Sec. V expresses our conclusions. 



II. BOUND STATE PION 

A. Bethe-Salpeter and gap equations 

Poincare covariance entails that the Bethe-Salpeter 
amplitude for an isovector pseudoscalar bound-state of 
a dressed-quark and -antiquark takes the form 

ri(fc; P) = T^75 [iE^{k; P) + 7 • P) 

+ -fk G^{k- P) + a^^k^PMk; P)] , (1) 

where k is the relative and P the total momentum of 
the constituents, and {T^,j — 1,2,3} are the Pauli ma- 
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trices.^ This amplitude is determined from the homoge- 
neous Bethe-Salpeter equation (BSE): 

[Ti{k-P)]tu = j ■0j^[xi{q;P)]srKi:{q,k;P) , (2) 

where xiiq; P) = Siq + P)ri{q;P)S{q), r,s,t,u rep- 
resent colour, flavour and spinor indices, and K is the 
quark-antiquark scattering kernel. In Eq. (2), S is the 
dressed-quark propagator; viz., the solution of the gap 
equation: 



s{py 



i'^ ■ p + m 



g Df,^{p - q)—-/^S{q)—T^{q,p), (3) 



wherein m is the Lagrangian current-quark mass, D^^, is 
the gluon propagator and Ti, is the quark-gluon vertex. 



B. Momentum-independent vector-boson exchange 



Following Ref. [22] we define 



(4) 



where tog is a gluon mass-scale (such a scale is gen- 
erated dynamically in QCD, with a value ^ 0.5 GeV 
[23]) and proceed by embedding this interaction in a 
rainbow-ladder truncation of the DSEs. This means 
^u{Pj q) — Iv in both Eq. (3) and the construction of K in 
Eq. (2). Rainbow-ladder is the leading-order in a nonper- 
turbative, symmetry-preserving truncation [24, 25]. It is 
known and understood to be an accurate truncation for 
pseudoscalar mesons [26, 27]. 

With this interaction the gap equation becomes 

S-\p) ^^^.p + n^+^-LJA^^^^ 5(g) 7, • (5) 



The integral possesses a quadratic divergence, even in the 
chiral limit. If the divergence is regularised in a Poincare 
covariant manner, then the solution is 



s{py 



i"f ■ p + M ^ 



(6) 



where M is momentum-independent and determined by 
M , 1 



M = 



M2 



(7) 



^ We employ a Euclidean metric with: {7p,7i/} = 25fj,y; =7^1; 
75 = 74717273; and a ■ b = 5I]i=i Qj^i- A timeliko four-vector, 
Q, has < 0. Furthemore, we consider the isospin-symmetric 
limit. 



To proceed we must specify a regularisation procedure. 
We write [28] 



S + M2 



2 

s + AP 



(8) 



(9) 



where uv are, respectively, infrared and ultraviolet reg- 
ulators. It is apparent from Eq. (9) that a nonzero value 
of Til- =: l/Aii- implements confinement by ensuring the 
absence of quark production thresholds [29, 30]. Further- 
more, since Eq. (4) does not define a renormalisable the- 
ory, Auv ^/t-on cannot be removed but instead plays 
a dynamical role and sets the scale of all dimensioned 
quantities. 

The gap equation can now be written 



M 



(10) 



where C'''{AP)/AP = r(-l,A/V2^) - r(-l,A/V2), 
with r(a, y) being the incomplete gamma-function. 

Using the interaction we've specified, the homogeneous 
BSE for the pseudoscalar meson is (q^ = q + P) 



r.(P) 



4 1 
3^ 



d^q 
(2^ 



7MX7r(g+,9)7M 



(11) 



With a symmetry-preserving regularisation of the inter- 
action in Eq. (4), the Bethe-Salpeter amplitude cannot 
depend on relative momentum. Hence Eq. (1) reduces to 



r.(P) = 75 



(12) 



Crucially, Fj^^P), a component of pseudovector origin, 
remains. It is an essential component of the pion, which 
has very significant measurable consequences and thus 
cannot be neglected. 



C. Ward-Takahashi identity 

Preserving the vector and axial- vector Ward-Takahashi 
identities is essential when computing properties of the 
pion. The to = axial- vector identity states 

P^n^{k+, k) = S-\k+)ij, + ^755^l(fc) ^ (-^3) 

where r5^(fc_)_, k) is the axial- vector vertex, which is de- 
termined by 



(14) 
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To achieve this, one must implement a regularisation that 
maintains Eq. (13). To see what this entails, contract 
Eq. (14) with and use Eq. (13) within the integrand. 
This yields the following two chiral limit identities: 



M 




8 M 



3 



1 



(2^ 



P ■ 



NP q\ 
P-q 



AP 



A/2 



AP 



(15) 
(16) 



which must be satisfied after regularisation. Analysing 
the integrands using a Feynman parametrisation, one ar- 
rives at the follow identities for P^ = = m: 



M 



16 M 



d-^q 



1 



AP 



(27r)4 [g2 + M2]2 



(17) 
(18) 



Equation (17) is just the chiral-limit gap equation. 
Hence it requires nothing new of the regularisation 
scheme. On the other hand, Eq. (18) states that the 
axial-vector Ward-Takahashi identity is satisfied if, and 
only if, the model is regularised so as to ensure there 
are no quadratic or logarithmic divergences. Unsurpris- 
ingly, these arc the just the circumstances under which a 
shift in integration variables is permitted, an operation 
required in order to prove Eq. (13). 

We observe in addition that Eq. (13) is valid for ar- 
bitrary P. In fact its corollary, Eq. (15), can be used 
to demonstrate that in the chiral limit the two-flavour 
scalar-meson rainbow-ladder truncation of the contact- 
interaction DSEs produces a bound-state with mass 
TOct ~ 2 M [31]. The second corollary, Eq. (16) entails 



da [C'"(cj(Af2,a,p2)) + C^\uj{AP ,a, P"^)] 



with uj{A/P,a,P^) = AP + a{l - a)P^ and Ci"(z) 
-z(d/dz)C™(z). 



(19) 



where 
ICee 



Kef — 
JCfe = 

JCpF = 



da [C'"(cj(A/2,a,-m^)) 



+2a{l~a)mlcT{io{AP,a,~ml)) , (21) 



-ml / daC^ {uj{AP,a,-ml)), 
Jq 

^APJ^ dat^{u{AP,a-ml)), 
—2Kfe 1 



(22) 



(23) 
(24) 



with Ci{z) = Ci{z)/ z. We used Eq. (19) to arrive at this 
form of JCff- 

In the computation of obscrvables, one must use the 
canonically-normalised Bethe-Salpeter amplitude; i.e., 
is rescaled so that 

= NctrJ ^T^{-P)-^S{q + P) T^{P) S{q) . 

(25) 

In the chiral limit, this means 

1 = ^^C^{A'P;tItI)E^[E^ - 2F^]. (26) 

With the parameter values (in GeV) [22] 

Tog = 0.11 , Air = 0.24 , Auv = 0.823 , (27) 

one obtains the results presented in Table I. We note 
that the Icptonic decay constant and in-pion condensate 
are given by the following expressions [32, 33]: 



M 2^2 JCIe 



-EE 



EF 



(28) 
(29) 



In the chiral limit ^ = ~{qq)\ i.e., the so-called 
vacuum quark condensate [34]. Moreover, also in this 
limit, one may readily verify that [22] 



m=0 



Ad 



(30) 



which is a particular case of one of the Goldberger- 
Treiman relations proved in Ref. [32]. 



D. Pion's Bethe-Salpeter kernel 

We are now in a position to write the explicit form of 
Eq.(ll): 



E^{P) 

FrriP) 



ICee Kef 
Kfe Kff 



E. [P) 

F. {P) 



(20) 



E. Dressed-photon-quark vertex 

In coupling photons to a bound-state constituted from 
dressed-quarks, it is important that the quark-photon 
vertex be dressed so that it satisfy the vector Ward- 
Takahashi identity [12]. Indeed, where possible it should 
be dressed at a level consistent with the truncation used 
to compute the bound-state's Bethc-Salpctcr amplitude 
[35] . With our treatment of the interaction described in 



TABLE I: Results calculated with the parameter values in 
Eq. (27): is obtained from Eq. (20); k^, /,r are defined in 
Eqs. (29), (28); nip is determined by solving Eq. (37); and the 
charge radii are discussed in connection with Eq. (38). The 
static properties are commensurate with results from QCD- 
based DSE studies [33]. (Dimensioned quantities are listed in 
GeV or fm, as appropriate.) 







M ^/l^ rup 


' q 'TV ' 7T 


m=0 
m=0.008 


4.28 0.69 
4.36 0.72 


0.40 0.22 0.094 0.90 
0.41 0.14 0.22 0.094 0.91 


0.34 0.30 0.45 
0.33 0.30 0.44 



connection with Eq. (4) . the bare vertex 7^ is sufficient 
to satisfy the Ward-Takahashi identity and ensure, e.g., 
a unit value for the charged pion's electromagnetic form 
factor [22]. However, given the simplicity of the DSE 
kernels, one can readily do better. 

A vertex dressed consistently with our rainbow-ladder 
pion is determined by the following inhomogeneous 
Bcthe-Salpeter equation: 



(27r)4 '"-^^ 



laXp.iQ+,<lha, (31) 



where X^^{q+,q) = Siq + P)T^{Q)S{q). Owing to the 
momentum-independent nature of the interaction kernel, 
the general form of the solution is 



r,iQ)=^lPTiQ') + 7^PLiQ^) 



(32) 



where Qf^jJ^ = and 7 J + 7^ = 7^- This simplicity 
doesn't survive with a more sophisticated interaction. 

Upon insertion of Eq. (32) into Eq. (31), one can readily 
obtain 

= (33) 
owing to Eq. (16). Using this same identity, one finds 

1 



1 + A%(Q2 



(34) 



with 



X / daail~a)Q^cT{uj{M^,a,Q'^)). (35) 



It is plain that 



(36) 



so that at Q2 = in the rainbow-ladder treatment of the 
interaction in Eq. (4) the dressed-quark-photon vertex is 
equal to the bare vertex. ^ However, this is not true for 
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FIG. 1: Dressing function for the transverse piece of the 
quark-photon vertex; viz., Pt{Q^) in Eq. (34). 



^ 0. Indeed, the transverse part of the dressed-quark- 
photon vertex will exhibit a pole at that < for which 



l + if-,(Q2) = 0. 



(37) 



This is just the model's Bethe-Salpeter equation for the 
ground-state vector meson. The mass obtained therefrom 
is listed in Table I. 

In Fig. 1 we depict the function that dresses the trans- 
verse part of the quark-photon vertex. The pole as- 
sociated with the ground-state vector meson is clear. 
Another important feature is the behaviour at large 
spacelike-Q^; namely, Pt{Q^) 1^ as — s- 00. This 
is the statement that a dressed-quark is pointlike to a 
large-Q^ probe. The same is true in QCD, up to the loga- 
rithmic corrections which are characteristic of an asymp- 
totically free theory [35]. 

One can define an electromagnetic radius for a dressed- 
quark; viz.. 



= -6 



(38) 



Our computed value is reported in Table I. It is notewor- 
thy that TqUipl = 0.63 so that, although the ground- 
state vector-meson pole is a dominant feature of Pt{Q^) 
in the vicinity of = 0, it does not completely deter- 
mine the electromagnetic radius of the dressed-quark. 

Nor, in fact, of anything else, as one can infer from the 
computed values of the pion charge radius reported in Ta- 
ble I:^ is obtained with r^((5) = 7/^; is computed 
with r^(Q) = 7m^t(Q'); and (r^)^ = + r^, which 
is just a consequence of the product rule. This empha- 
sises again that single-pole vector-meson-dominance is a 
helpful phenomenology but not a hard truth [5, 6, 35]. 



^ Equations (33), (36) guarantee a massless photon and demon- 
strate that our regularisation also ensures preservation of the 
Ward-Takahashi identity for the photon vacuum polarisation 
[36]. 



^ In computing Ttt , we follow Rcf. [22] . Owing to the vector Ward- 
Takahashi identity, the longitudinal part of the vertex does not 
contribute to the pion's elastic form factor. 



5 



3 - 



E 



p - pole arising from 

dressing tine quark-photon vertex 



4 6 

;2\ 



10 



FIG. 2: -F™(Q^) computed in rainbow-ladder truncation 
from the interaction in Eq. (4): solid curve - fully consistent, 
i.e., with a dressed-quark- photon vertex so that the p-pole 
appears; and dashed curve - computed using a bare quark- 
photon vertex. Dotted curve ~ fit to the result in Ref. [37], 
which also included a consistently-dressed quark-photon ver- 
tex and serves to illustrate the trend of contemporary data. 



We show in Fig. 2 that dressing the quark-photon 
vertex does not qualitatively alter the behaviour of 
at spacelike momenta. In particular, it does 
not change the fact that a momentum-independent in- 
teraction, Eq. (4), regularised in a symmetry-preserving 
manner, produces^ 



A. Anomaly 

We first consider the chiral limit and ~ 0, in 
which case Eq. (42) describes the "triangle diagram" that 
produces the Abelian anomaly and one must compute 
G(0,0, 0). We have explained above that our regulari- 
sation of the interaction in Eq. (4) ensures that the non- 
anomalous vector and axial- vector Ward-Takahashi iden- 
tities are satisfied. The outcome for the anomalous case 
is therefore very interesting. 

Two contributions are obtained upon inserting Eq. (12) 
into Eq.(42); viz., one associated with E.^{P), which we'll 
denominate Ge, and the other with FTr{P), to be called 
Gp- We first examine the latter. To obtain Gi?(0, 0, 0) 
one need only expand the integrand in Eq. (42) around 
ki = = k2 and keep the term linear in F.^(P)kiak2p, a 
process which yields 

Gf(0,0,0) = -^ / dss''F^{P)av{s) 
X [av(s)2 -I- sav{s)a'y{s) + crs(s)cr^(s)] , (44) 
where (t'{s) = ^"'(•s) and we have written 

Sii) = -i-fiavie^) + <Js{e^). (45) 

Using Eq. (6), one readily finds that a'y = —Cy, a'g = 
—May. These identities, when inserted into Eq. (44), 
reveal that the integrand is identically zero, so that 



FT\Q^ ^ oo) = constant. 



(39) 



Gf(0,0,0) = 0. 



(46) 



III. TRANSITION FORM FACTOR: 7*77^7 

In the rainbow-ladder truncation this process is com- 
puted from [16] 

%y{ki,k2) =Tf_,^iki,k2) + T^f^ik2,ki) , (40) 

where the pion's momentum P = ki + k2, ki and fc2 are 
the photon momenta and 

T)ii/(fcl,fc2) = e.fj,yal3klak2l3G{kl,ki -^2,^1) (41) 

(2^ 



Wtt^z^i xAh,i2) «Qr^(^2, ^12) s{h2) iQT Aei2,ii), 



(42) 

with £1 = £ ~ ki, £2 = £ + k2, li2 = £ - ki + k2, and 
Q = diag[e„,ed] = ediag[2/3, -1/3], aom = e2/(47r). 
The kinematic constraints are: 

fc? = Q2 fc| = , 2 fci • fc2 = -(m^ + q2) (43) 



* The rainbow-ladder truncation omits a so-called meson-cloud 
component of F^™(Q^) but this, too, affects only the behaviour 
in a measurable neig hbourhood of = [38] . 



This is a particular case of the general result proved 
in Ref. [14]. As explained therein, since the integral in 
Eq. (44) is logarithmically divergent, the result is only 
transparent with the choice of momentum-partitioning 
that we have employed. 

The remaining contribution is G^;, which, following the 
methods of Sec. II, can be written 



G£(0,0,0) 



dHE^{P)av{£l2)ay{£l)av{tl). 



(47) 

The integral is convergent and therefore a shift in inte- 
gration variables cannot affect the result. It follows that 



G£(0,0,0) = £;,(P)^ / ds:. 



(S + M2)3 



(48) 



If we employ Eq. (8), as with all other computations 
hitherto, this becomes 



GB(0,0,0) = ^Ci"(M2), 



(49) 



where 4"(z) = {z^ /2){d'^ / dz^)C'''{z) and we have used 
the Goldberger-Treiman relation in Eq. (30). 

In what has long been a textbook result, the anoma- 
lous Ward-Takahashi identity states that Ge{^, 0, 0) = ^: 
truly, just this simple fraction. Equation (49) is plainly 
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inconsistent with this because it produces a number that 
depends on the values of the parameters Ajr, Auv- Indeed, 
with the values in Eq. (27), our regularisation of Eq. (4) 
gives /^G£;(0,0,0) = 0.36. What has gone wrong? 
The answer lies in the observation that 

°(A/2) ^ C2{M\ tI ^ oo, tI, ^ 0) = i Af2 . (50) 

One could have judged at the outset that no regulari- 
sation scheme which bounds the loop-integral can sup- 
ply the correct result for the anomalous Ward-Takahashi 
identity because it blocks the crucial connection between 
the anomaly, topology and DCSB [39]. 

To elucidate, return to Eq. (48). The integral is conver- 
gent and dimensionless. Hence, it cannot depend on M . 
In a particular application of the procedure elucidated 
in Refs. [10, 12-14], the change of variables C(s) = M/s 
yields 

G£;(0,0,0) = ^ dCj^^-^ (51) 
1 f°° d 1 1 

The last line emphasises the connection between the sim- 
ple rational-number result and the spacetime boundary: 
the anomaly is determined by the integral of a total 
derivative. The result in Eq. (52) is obtained if, and only 
if, chiral symmetry is dynamically broken, since in this in- 
stance alone can Eq. (30) be used to completely eliminate 
the pion structure- factor: Ej^iP), from the expression. 

B. Asymptotic behaviour 

1. Massless vector-boson exchange 

In Ref. [40] , using the methods of light- front quantum 
field theory, it was shown that 

hm Q^G{Q\^\Q\{)) = A^^fl (53) 

It is notable that this is a factor of Tr/as{Q'^) bigger than 
the kindred limit of the elastic pion form factor [40-42]; 
i.e., at = 4GeV^, more than an ordcr-of-magnitude 
larger. 

Our analysis of Eq. (42) is performed within a Poincarc 
covariant formulation. In this elucidated in 

Ref. [20] , the asymptotic limit of the doubly-off-shell pro- 
cess (7*7* — ?> tt") is reliably computable in the rainbow- 
ladder truncation because both quark-legs in the dresscd- 
quark-photon are sampled at the large momentum-scale 
Q^, with the resuh [18, 19] 

lim Q'GiQ^-Q'-ml/2,Q')^li7T'fl (54) 

if the propagator of the exchanged vector-boson behaves 
as 1/fc^ for large-A:^. In order to obtain this result it 



is crucial that the pion's Bethe-Salpeter amplitude de- 
pends on the magnitude of the relative momentum and 
behaves as for large-fc^, as it does in QCD. (See also 
Sec. inc.) 

Equations (53) and (54) correspond to the asymptotic 
limits of different but related processes. Part of the 
mismatch owes to the fact that in the process 7*77r° 
not all quark-legs attached to vertices carry the large 
momentum-scale Q^; namely, £2 in Eq. (42), and hence 
some amount of vertex dressing contributes, even at large 
Q^. This is consistent with a more general observa- 
tion; namely, that in a covariant calculation any number 
of loops contribute to 7*77r° at leading order [40], and 
these provide a series of logarithmic corrections which 
should be summed. The same is true of the pion's elas- 
tic form factor [14]. Nevertheless, the correct power-law 
behaviour is necessarily produced. 



2. Contact interaction 

With the QCD-based expectation made clear, we now 
turn to the outcome produced by the contact interac- 
tion, Eq. (4). In this instance the arguments used to ob- 
tain Eq. (54) fail conspicuously because the pion's Bethe- 
Salpeter amplitude is completely independent of the rel- 
ative momentum: all values of the relative momenta are 
equally likely. This is why the interaction yields Eq. (39), 
a result in striking disagreement with experiment. A sim- 
ilar result is obtained in the present context. However, a 
decision must be made before that can be exhibited. 

Recall Eqs. (44) and (48): the first is logarithmically 
divergent while the second is convergent even if the reg- 
ularisation parameters are removed. Indeed, one needs 
to remove the regularisation scales if the anomaly value 
is to be recovered. However, the form factor is then ill- 
defined because the F7r(P)-term contributes the logarith- 
mic divergence just noted. We proceed by removing the 
regularisation in computing Ge but retaining it in cal- 
culating Gf- Notably, as we will see, with Eq. (4) no 
internally consistent scheme can provide QCD-like ultra- 
violet behaviour but this prescription serves to preserve 
the infrared behaviour. 

There is one more step in implementing this scheme. 
In arriving at expressions such as those defining the 
pion's Bethe-Salpeter kernel (see Sec. II D), we re-express 
a product of propagator-denominators via a Feynman 
parametrisation, then perform a change-of-variables, and 
thereafter rewrite the result using Eq. (8). This does 
not introduce any difficulties when the boundary at 
spacctime-infinity has no physical impact. However, as 
we have seen, that is not the case with the anomaly. The 
integral which defines Gf is logarithmically divergent. A 
shift of integration variables changes its value, and in do- 
ing that one runs afoul of the fact that it is impossible 
to simultaneously preserve the vector and axial vector 
Ward-Takahashi identities for triangle diagrams in field 
theories with axial currents that are bilinear in fermion 
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FIG. 3: (Colour online) 7*7r"7 transition form factor in 
Eq. (42). Solid curve - Full computation described in text, 
including -Fir(-P) contribution; dashed curve - result obtained 
without dressing the quark-photon vertex; dash-dotted curve 
- result obtained with Ft^IP) = 0, i.e. forced artificially to 
vanish; dotted curve - monopole fit to the QCD-based result 
in Ref. [20], which agrees with the data reported in Refs. [2, 3], 
green diamonds and blue squares, respectively. 



fields. Any shift in variables from that used in Eq. (42) 
changes the value of 6*^(0,0,0). We compensate by an 
additional regularising subtraction; i.e., by redefining 



> 

CD 




-{ml 
\ml 



)2)/2,0)-Gf(0,0,0). (55) 
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FIG. 4: (Colour online) 7*7 — >■ transition form fac- 
tor. Data: red circles, Ref. [1]; green diamonds, Ref. [2]; 
and blue squares, Ref. [3]. Solid curve - G((5^,0) computed 
using the symmetry-preserving, fully-self-consistent rainbow- 
ladder treatment of the contact interaction in Eq. (4), with the 
dotted-curve at right showing its nonzero asymptotic limit; 
and dot-dashed curve ~ fit to the 7*7 — >■ n° transition form 
factor computed in a QCD-based rainbow-ladder-truncation 
DSE study [20]. The curves have been divided by f-^) in 
order to match the normalisation of the data. 



In doing so we implement an anomaly-free electromag- 
netic current [43]. 

In Fig. 3 we depict the result produced from Eq. (4) 
using the regularisations just described. Comparing the 
solid- and dashed-curves, it is evident that the effect of 
dressing the quark-gluon vertex diminishes with increas- 
ing and therefore has no impact on the asymptotic 
behaviour of the transition form factor. It does, however, 
affect the ncutral-pion interaction-radius, which can be 
defined via 



&^\nG{Q\-{ml + Q^)l2,Q) 



(56) 



This yields r*o ~ 0.30 fm and rja = 0.45 fm, values that 
are not sensibly distinguishable from the chargcd-pion 
values listed in Table I. Near equality of and r*o 
is also found in the QCD-based calculations of Refs. [20, 
35]. 5 

More significantly, the solid- and dashed-curves in 
Fig. 3 show that, as with the elastic form factor [22], the 



^ Choosing instead the 7*7* — > tt" form factor, one finds r*o = 
0.44 fm and r^g'' = 0.55 fm, values which are larger because the 
momentum-scale enters into both quark-photon vertices. 



presence of the pion's necessarily-nonzero pseudovector 
component, F7r(P), leads to 



lim G(Q^ 



-{ml + Q^)/2, 0) = constant. (57) 



This is consistent with the picture developed in Ref. [22]; 
namely, it is possible to treat the contact interaction, 
Eq. (4), so that it yields static properties of the pion in 
agreement with experiment and computations based on 
well-defined and systematically improvable truncations 
of QCD's DSEs. However, a marked deviation from ex- 
periment occurs in processes that probe the pion with 
> and it is impossible to obtain results which 
agree with perturbative-QCD, even at the gross level of 
form-factor power-laws. 

These observations are emphasised by the comparisons 
presented in Figs. 4 and 5. The 7*7 — > tt" form fac- 
tor obtained using the symmetry-preserving, fully-self- 
consistent rainbow-ladder treatment of the contact in- 
teraction in Eq. (4) is in glaring disagreement with all 
existing data. This is what it means to have a pointlike- 
component in the pion: all form factors must asymp- 
totically approach a constant. That limit rapidly be- 
comes apparent with increasing momentum transfer be- 
cause the dynamically generated mass-scale associated 
with low-energy hadron phenomena is M ^ 0.4 GeV. No 
study that neglects the pion's pseudovector component 
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0.10 



0.05 



0.00 



MX) = constant 




';,(X)~X(1-X) 



It follows from this expression that 
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Cf(GeV2) 
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FIG. 5: (Colour online) Q^-weighted 7*7 — >■ tt" transi- 
tion form factor. Data: red circles, Ref. [1]; green dia- 
monds, Ref. [2]; and blue squares, Ref. [3]. Solid curve - 
Q^G{Q^ ,0) computed using the symmetry-preserving, fully- 
self-consistent rainbow-ladder treatment of the contact in- 
teraction in Eq. (4), which produces = constant (see 
Sec. Ill C; and dot-dashed curve - fit to the 7*7 — ^ tt" tran- 
sition form factor computed in a QCD-based rainbow-ladder- 
truncation DSE study [20]. Both curves have been divided by 
(27r^/^) in order to match the data's normalisation. 



can provide a valid explanation or interpretation of the 
7*7 7r° transition form factor, or any other of the 
pion's form factors. 

On the other hand, it is noteworthy that the DSE re- 
sult [20] , which is based on an interaction that preserves 
the one-loop renormalisation group behavior of QCD, 
agrees with all but the large-Q^ BaBar data. 



C. Pion distribution amplitude 

It is worthwhile to consider a little further the nature 
of a pointlike pion. As explained in Ref. [44], with the 
dressed-quark propagator and pion Bcthe-Salpeter am- 
plitude in hand, one can compute the pion's valence- 
quark parton distribution function in rainbow-ladder 
truncation. For the contact interaction, the result is 

3 f d^i 

X r^{-P)S{£)n-^S{£)r^{P)S{£~P), (58) 

where v? = Q, n - P — P+, and x is the Bjorken variable. 



{n-P) 



n+l 



2i 



tro 



-(n-£)" 



X r^{-P)Sie)n-jSi£)T^{P)S{£-P). (59) 

At this point we'll specialise to the chiral limit and: eval- 
uate the Dirac-trace; use a Feynman parametrisation to 
re-express the product av{£'^)crv{{£ — P)^) which arises; 
shift variables £ ^ {£ + aP), where a is the Feynman 
parameter; use the 0(4) invariance of the measure to 
evaluate the angular integrals; and thereby arrive at 



/ dxx"qv{x) ~ — 
Jo " 



1 



n-f 1 ' 



(60) 



where the last line follows because the pion's Bethe- 
Salpeter amplitude is canonically normalised, Eq. (26). 

The distribution function is readily reconstructed from 
Eq.(60); and one finds that even with inclusion of the 
pion's necessarily- nonzero pseudovector component, the 
contact-interaction produces 

qvix) = e{x)e{i - x) , (61) 

which corresponds to a pion distribution amplitude 



(f).^{x) 



constant. 



(62) 



This outcome provides another way of understanding the 
inability of the contact interaction to reproduce the re- 
sults of QCD. 

As reviewed and explained in Ref. [44], Goldstone's 
theorem in QCD is expressed in a remarkable corre- 
spondence between the quark-propagator and the pion's 
Bethe-Salpeter amplitude; i.e., between the one- and 
two-body problems [32]. The long- known fact that the 
dressed-quark mass function behaves as [45-49] 



large-p^ 1 
p2 



(63) 



entails that in QCD every scalar function in the pion's 
Bethe-Salpeter amplitude, Eq. (1), depends on the rela- 
tive momentum, fc, as ~ for large-fc^ (with addi- 
tional logarithmic suppression). It is impossible to find 
a kinematic arrangement of the dressed-quarks consti- 
tuting the pion in which the Bethe-Salpeter amplitude 
remains nonzero in the limit of infinite relative momen- 
tum. 

It follows that in QCD the pion's valence-quark distri- 
bution behaves as (1 — x)'^'^'^ , < 7 ^ 1, for a; ~ 1, at a 
renormalisation scale of 1 GeV [44, 50] . Hence there is no 
renormalisation scale in the application of perturbative 
QCD at which Eq. (61) is a valid representation of non- 
perturbative QCD dynamics; namely, no scale at which 
it is tenable to employ (t>Ti{x) = constant, or even, more 
weakly, fiat and nonzero at x = 0, 1. 



9 




FIG. 6: (Colour online) Q^-weighted 7*7 — >■ tt" transition 
form factor. Data: red circles, Ref. [1]; green diamonds, 
Ref. [2]; and blue squares, Ref. [3] - all normalised by the 
asymptotic form in Eq. (53) so that the pQCD limit is marked 
by the dotted line at "1". Solid curve - Q^Gi^(Q^O) cal- 
culated as described herein, normalised by the asymptotic 
form in Eq. (64); dashed curve - monopole with mass-scale 
(2/3)47r'^/^, which bounds from above the 7*7* — >■ vr" tran- 
sition form factor computed in Ref. [20], normalised by the 
asymptotic form in Eq. (54); and dot-dashed curve - monopole 
with mass-scale (7/8)47r^/.^, obtained as fit to the 7*7 — >■ tt" 
form factor computed in Ref. [20], normalised by this same 
mass-scale. 



tion yields 
G^(Q^-(g2-fm^)/2,0) Q'>1'''^_^ 



In 



, (64) 



cf. Eq. (54). A similar result is obtained for the doubly- 
ofF-shell process, with the only difFcrcncc being that the 
power of the logarithm is reduced "2 — > 1" . 

We stress that in Eq. (64), AI is the dressed-quark 
mass. Table I emphasises that M is a computed quan- 
tity, which is completely determined once the interac- 
tion and truncation are specified. The value of M is 
tightly connected with those of all other measurable 
properties in the Table. Thus, in a well-constrained and 
internally-consistent analysis, one cannot significantly al- 
ter M without materially changing all the other static 
properties which characterise the pion. No theoretical 
analysis is reliable if it allows itself to skirt these con- 
straints. 

In Fig. 6 we depict the contact-interaction result for 
Q^Ge{Q^ ,~{Q^ + toJ)/2,0) normalised by the asymp- 
totic form in Eq. (64). In addition, we plot a monopole 
with mass-scale (2/3)47r^/^, which bounds uniformly 
from above the QCD-based DSE calculation of the 
7*7* — >■ tt" transition form factor reported in Ref. [20]; 
and a monopole with mass-scale (7/8)47r^/^, which is a 
fit to the 7*7 — >■ TT*' transition form factor also computed 
therein. The origin of these mass-scales was discussed in 
Sec. Ill Bl. It is striking that these curves all approach 
their asymptotic limits from below. Stated differently, 
each is a monotonically-increasing concave function. In- 
deed, this is even true of the solid curve in Fig. 5. 



IV. REFLECTIONS ON EXTANT DATA 

We have shown that an internally consistent treat- 
ment of the contact interaction is incompatible with ex- 
tant pion elastic (Fig. 2) and transition form factor data 
(Figs. 4, 5). Notwithstanding this, the results elucidated 
can be used in combination with QCD-based DSE stud- 
ies in order to comment on available data for the process 
7*7 -J> tt". 

To begin we remark upon a similarity between the Q^- 
dependcncc of the dash-dot- and dottcd-curves in Fig. 3; 
i.e., the QCD-based DSE result and that obtained from 
the contact interaction if the pion's pseudovector com- 
ponent is artificially eliminated. We emphasise that if 
Ft^^P) is forced to zero, then one is no longer represent- 
ing faithfully the features and consequences of a vector- 
vector contact interaction. Hitherto, this has neverthe- 
less been a conventional mistreatment of Eq. (4). Its 
consequences were first elucidated in Ref. [22]. We de- 
scribe results obtained through this intervention in order 
to clarify its real implications for the 7*7 — > tt" transition 
form factor. 

An analysis of Eq. (42) shows that the contact interac- 



V. CONCLUSIONS 

We have shown that in fully-self-consistent treat- 
ments of pion: static properties; and elastic and tran- 
sition form factors, the asymptotic limit of the prod- 
uct Q^G-y.-yTrO (Q^), which is determined a priori by 
the interaction employed, is not exceeded at any finite 
value of spacelike momentum transfer: the product is a 
monotonically-increasing concave function. We under- 
stand a consistent approach to be one in which: a given 
quark-quark scattering kernel is specified and solved in 
a well-defined, symmetry-preserving truncation scheme; 
the interaction's paramcter(s) are fixed by requiring a 
uniformly good description of the pion's static proper- 
ties; and relationships between computed quantities are 
faithfully maintained. 

Within such an approach it is nevertheless possi- 
ble for Q2^™(Q2)^ with i^™(Q2) being the elastic 
form factor, to exceed its asymptotic limit because 
the leading-order matrix-element involves two Bethe- 
Salpeter amplitudes. This permits an interference be- 
tween dynamically-generated infrared mass-scales in the 
computation. Moreover, for F™((3^) the perturbative 
QCD limit is more than an order-of-magnitudc smaller 
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than m^. Owing to the proximity of the p- meson pole 
to (5^ = 0, the latter mass-scale must provide a fair 
first-estimate for the small-Q^ evolution of F°™{Q'^). A 
monopole based on this mass-scale exceeds the pQCD 
limit VQ^ > 0. For the transition form factor, however, 
the opposite is true because is less-than the pQCD 
limit, Eq. (53). 

A vector current-current contact-interaction may be 
described as a vector-boson exchange theory with vector- 
field propagator (1/fc^)'^, k = 0. We have shown (see 
Figs. 4. 5) that the consistent treatment of such an inter- 
action produces a 7*7 — > tt" transition form factor that 
disagrees with all available data. On the other hand, 
precisely the same treatment of an interaction which pre- 
serves the one-loop renormalisation group behaviour of 
QCD, produces a form factor in good agreement with all 
but the large-Q^ data from the BaBar Collaboration [1]. 

Studies exist which interpret the BaBar data as an 
indication that the pion's distribution amplitude, 4>.^{x), 
deviates dramatically from its QCD asymptotic form, in- 
deed, that 4>-n{x) = constant, or is at least flat and nonva- 
nishing at a; = 0, 1 [51, 52]. We have explained that such 
a distribution amplitude characterises an essentially- 
pointlike pion; and shown that, when used in a fully- 
consistent treatment, it produces results for pion elas- 
tic and transition form factors that are in striking dis- 
agreement with experiment. A bound-state pion with 
a pointlike component will produce the hardest possible 
form factors; i.e., form factors which become constant at 
large-Q^. 

On the other hand, QCD-based studies produce soft 



pions, a valence-quark distribution amplitude for the pion 
that vanishes as (1 — x)^ for x ^ 1, and results that 
agree well with the bulk of existing data. 

Our analysis shows that the large-Q^ BaBar data is in- 
consistent with QCD and also inconsistent with a vector 
current-current contact interaction. It supports a conclu- 
sion that the large-Q^ data reported by BaBar is not a 
true representation of the 7*7 — >■ tt*^ transition form fac- 
tor, a perspective also developed elsewhere [53]. We are 
confirmed in this view by the fact that the 7* — 777 and 
7* — > rj'^ transition form factors have also been measured 
by the BaBar Collaboration [54], at 112GeV^ and 

in these cases the results from CLEO [3] and BaBar are 
fully consistent with perturbative-QCD expectations. 
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